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We apply a diagrammatic expansion method around the atomic limit (U >> t) for the U-t-t' 
Hubbard model at half filling and finite temperature by means of a continued fraction representation 
of the one-particle Green's function. From the analysis of the spectral function A(k, uf) we find an 
energy dispersion relation with a (cos/d — cosfcy) modulation of the energy gap in the insulating 
phase. This anisotropy is compared with experimental ARPES results on insulating cuprates. 
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During the last decade of research work on high tem- 
perature superconductors angular resolved photoemis- 
sion spectroscopy (ARPES) has played an important role 
in elucidating their electronic properties An under- 
standing of the single particle properties in the normal 
state is a prerequisite for a theory of a mechanism for su- 
perconductivity as well as for transport properties. Over 
the years ARPES data have continuously provided sur- 
prises and new insights and thereby served as a guid- 
ance to theoretical developments. The observation of the 
cL^^ -shape of the energy gap, pseudogap structures in 
the metallic phase [|]||], strong anisotropies in the quasi- 
particle (qp) peak lineshapes and a possible partial de- 
struction of the Fermi surface 0], or the unusual fre- 
quency and temperature dependence of the qp peakwidth 

fare examples for intriguing information obtained from 
RPES experiments. 

For the normal state an important goal is the descrip- 
tion for the evolution with doping from the antiferromag- 
netic (AF) and insulating parent compounds to the over- 
doped superconductors. This demands control over the 
spectral features of the Mott-Hubbard insulating state as 
a starting point. Yet, it proved to be difficult to repro- 
duce ARPES data for the single hole dispersion in the 
insulating cuprate Sr2Cu02Cl2 Q within t-J or Hub- 
bard models. In particular, for momenta along the Bril- 
louin zone (BZ) axis next-nearest neighbor (nnn) or even 
longer range hopping amplitudes had to be introduced to 
achieve a reasonable comparison to the measured spec- 
tra ]7|"[l0||. Even more striking are recent results on AF 
Ca2Cu02Cl2 that demonstrated an anisotropy of the in- 
sulating energy gap which was claimed to follow closely a 
e£j.2_ a 2-wave modulation along a remnant Fermi surface 
with a modulation amplitude of ~ 300meV comparable in 
magnitude to the AF exchange interaction Jll],[l2| . It has 
been pointed out that a d x 2_ y 2 gap modulation might fol- 
low naturally in the context of a projected SO(5) theory 
p3[ unifying antiferromagnetism and d-wave supercon- 
ductivity by a symmetry principle jl4| . 

In this paper we show that a 7^ (k) = (cos k x — cos k y ) 2 
modulation of the AF energy gap is realized in the half- 
filled Hubbard model on a square lattice with nearest 



(nn) and nnn hopping amplitudes. This result is ob- 
tained in an analytic strong coupling expansion around 
the atomic limit following a recently proposed strategy 
with a mapping to a Jacobi continued fraction represen- 
tation for the propagator Given that the Mott- 
Hubbard insulator is the appropriate starting point for 
studying the hole doping evolution in cuprates this intrin- 
sic anisotropy may bear an important preformed struc- 
ture for the doped metallic phase. 

The single particle Matsubara Green's function for the 
Hubbard model is represented in terms of Grassmann 
fields 7* and 7 in the Feynman path integral representa- 
tion by 



G ij (ra\r'a') = -(T T c i(T (T)cUr')) 



= ~z j d ^ llcrT 1 * ja ' T ' exp (" S '[ 7 *' 7 ^ 

with the partition function Z and the action S = Ski 
Satom where 
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Satom [7* , 7] = J dT (q^ - A*) 



(3) 



Here, cJ CT creates an electron at site i with spin er, and we 
restrict the hopping amplitudes Uj to nn (t) and nnn (t') 
sites. U is the on-site Coulomb repulsion, fi the chemical 
potential, and f3 = 1/T the inverse temperature. 

Gij{Ta\r' a 1 ) is evaluated diagrammatically in terms of 
a cumulant expansion around the atomic limit, i.e. an ex- 
pansion in powers of t/U ifTq] . As early on suggested by 
Sarker p7| it is convenient to introduce auxiliary Grass- 
mann fields {ip* aT , ipia T } by performing a Hubbard-Stra- 
tonovich transformation with respect to the kinetic en- 
ergy term Eq. 0. From the Green's function for the 
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FIG. 1. Irreducible diagrams contributing to G(k,iui n ) of 
the U-t-t' Hubbard model up to the order (t/U) 3 . The ver- 
tices (•) represent local n-particle cumulants Gq . Each of 
the lines represents a hopping process with amplitude either 
t or t'. 



auxiliary fields Vij(T<r\i J a') = —(T T ip i(JT ip* a , T ,) the prop- 
agator Gij \{t<j\t' a') in Eq. Ill is obtained via the relation 
Gij(Ta\r'a') = (T' 1 ~ V) ^(tct|tV) g§. Here T de- 
notes the self energy and V is the non- interacting Green's 
function for the auxiliary fields, i.e. V = V + VTV. 

As is known from previous work on strong-coupling 
expansions for the Hubbard model [ pii||Ts| ] the Green's 
function calculated this way does not have the correct 
analytic properties, because higher order cumulants con- 
tain also higher order poles. This problem was circum- 
vented in Ref. pl| by mapping the Green's function as 
calculated by the diagrammatic expansion to a Green's 
function Gj in a finite Jacobi continued fraction repre- 
sentation [E3i 
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such that Gj has the same series expansion as G to the 
same order in t/U. 

Specifically, we have calculated the Fourier trans- 
formed Green's function G(k,iui n ) at half filling /i = ^ 
up to order (t/U) 4 . While our independently obtained 
results agree with Ref. fTH for nn hopping only, in the 
case of a finite nnn hopping t' an additional tt' 2 /U 3 dia- 
gram contributes to G(k, ito n ) as shown in Fig. |l|. In this 
type of diagrams the vertices represent local n-particle 
cumulants Gq (connected Green's functions). Each line 
between two vertices represents a hopping process be- 
tween two sites. In the absence of a Wick theorem for the 
local Green's functions the higher order cumulants must 
be calculated separately. The resulting algebraic expres- 
sions become very involved for higher order cumulants; 
for the performance of this diagrammatic expansion we 
have therefore developed a special purpose computer al- 
gebra code. 
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FIG. 2. Spectral function A(k, w) of the half-filled U-t-t' 
Hubbard model as a function of u along the path k = (0, 0) 
— > (n, 7r) — * (tt, 0) — » (0, 0) through the first Brillouin zone 
for U = Wt, T = 0.2t and t' = -0A5t as obtained from the 
strong coupling expansion to order (t/U) 4 . 



As a result we obtain a Jacobi continued fraction ex- 
pression for Gj that has eight fraction levels, i.e. the 
coefficients ajv in Eq. || vanish for N > 8. This termi- 
nation for Gj translates into eight simple poles. Each 
of the continued fraction coefficients a,, 6, (i = 1, . . . , 8) 
is given by a fourth order polynomial in t/U or t'/U, 
respectively, with coefficients that depend on T, fi, U, 
and k. The explicit result for Gj is accessible electroni- 
cally |19|. In Fig. we show the corresponding spectral 
function 
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for different k = (k x , k y ) along the path (0, 0) — > (it, tt) — > 
(tt, 0) — > (0,0) in the first BZ. On the scale of this plot 
only six out of eight poles of Gj(k,co) carry significant 
and visible spectral weight. 

In analyzing the spectrum we map out the dispersion 
of the lowest energy hole excitation, i.e. the dispersion 
of the first peak below the gap (marked by an arrow in 
Fig. ^|). Its spectral weight is largest at the BZ center 
and drops monotonically towards (tt, tt). Fig. || shows 
the peak dispersion for different temperatures. Upon 
cooling the system approaches the AF ordered ground 
state with a doubled unit cell and a reduced magnetic 
BZ (determined by cosfc^ + cosk y > 0). The dispersion 
along the zone diagonal approaches a perfectly symmetric 
shape with respect to the point kd = (§, \ ) on the mag- 
netic BZ boundary reflecting the growing AF spin cor- 
relations. Remarkably, along the BZ axis the dispersion 
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FIG. 3. Temperature dependence of the dispersion E(k) of 
the low energy peak (marked by an arrow in Fig. ^) along 
a selected path in the first Brillouin zone for U = lOt and 
t' = -0.45t. 

remains very fiat at all temperatures. In fact, the disper- 
sion along the axis is obtained flatter than in previous 
numerical studies of Hubbard or t-J models ]7|j9|, |20| , ^lf . 
In units of the exchange coupling J = At 2 /U the total 
bandwidth of the energy dispersion is 1.58 for T = Q.2t, 
which is roughly consistent with experiment and previous 
calculations. 

In Fig. | we compare our result for the single hole dis- 
persion with the ARPES data on Sr2Cu02Cl2 from Ref. 
P|. For U = lOt the parameters t and t' were chosen to 
obtain a best fit to the data. The overall agreement, in 
particular along the BZ axis, is quite satisfactory. Hop- 
ping amplitudes beyond nnn hopping are found unnec- 
essary for reproducing the flat dispersion along the BZ 
axis. The deviation of the theoretical result from the ex- 
perimental data points along the BZ diagonal, i.e. the 
lack of symmetry with respect to - as realized in the 
AF state - is due to the finite temperature. In fact, the 
range of applicability of the t/U expansion has a lower 
bound in temperature. (Following the arguments of Ref. 
|15| we estimate that our results are valid for T > 0.16t.) 
When the magnetic correlation length is much larger 
than the hopping range as given by the order of the t/U 
expansion, an accurate description of the single particle 
Green's function is no longer expected. The data were 
taken 100K above the Neel temperature Tjv = 256K of 
S^CuO^CFj; at this temperature the magnetic correla- 
tion length is already as large as 250A as measured by 
neutron scattering f22|j . Nevertheless, the flatness of the 
dispersion along the BZ axis in the t/U expansion is a 
robust and temperature insensitive feature. 



FIG. 4. A comparison with the ARPES dispersion for 
Sr2Cu02Cl2 taken from Ref. [6] (black diamonds) and Ref. 
[12] (white circles). Error bars are shown at some selected 
experimental data points. The parameters in the theory were 
chosen as U = 10*, T = 0.2t, t' = -0.45i, and t = 0.605eV. 

In Fig. U the calculated energy dispersion E(kd) — 
E(k) of the lowest energy peak and the experimental data 
are plotted along the BZ path from ka to (tt, 0) which is 
parametrized in the form 7d(fc)|/2. Originally, the data 
of Ref. |ll| (black squares) were anticipated to imply a 
d x 2 _ y 2 -modulation of the energy gap which in the para- 
metrization of Fig. ^ would translate into a straight line. 
The perfectly linear relation between E(kd) — E(k) and 
7d(fc)| along the chosen path would furthermore imply 
an unreasonable cusp- like feature of the dispersion on the 
magnetic BZ boundary. Indeed, the more recent data 
by Ronning et al. Jl2| (white circles in Fig. ||) resolve 
this problem and rather provide evidence for a quadratic 
dependence of the energy gap on |7d(fc) in the vicinity of 
fcrf. Given the experimental error bars of Ref. |(| included 
in Fig. |B| and in the absence of an error estimate for the 
recent data by Ronning et al. JX2J our results are clearly 
compatible with experiment. With the same parameter 
set used for the fit of the ARPES dispersion in Fig. || 
the gap modulation amplitude ~ 300meV as measured 
in Ca 2 Cu02Cl2 [|llj is reproduced as well. 

Fig. |6| shows the temperature dependence of the gap 
modulation. In fact, when plotted versus 7^(fc) it be- 
comes evident that the gap approaches a perfect 7^(fc) 
momentum dependence at low temperatures. For t' = 
the modulation vanishes in the expansion up to the order 
(t/U) 4 because the k dependence of all diagrams to this 
order arises in the form cos k x + cos k y for t' = which 
vanishes on the BZ boundary for t' = 0. If diagrams of 
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FIG. 5. Same data as in Fig. ^ along the BZ path from 
k d = (71-/2,71-/2) to (tt,0), but here E(k d ) - E(k) is plot- 
ted as a function of |7d(fc)|/2. In addition, ARPES data on 
Ca2Cu02Ci2 (from Ref. [11]) are shown (black squares). 



FIG. 6. Dispersion along the path from kd = (7r/2,7r/2) 
to (7T, 0) relative to its value at kd as a function of 
(cos fen, — cosfc a ) 2 /4 for different temperatures with U = Wt 
and t' = -OAt. 



order (t/U) 6 or higher were taken into account, a small 
modulation is expected to appear also for t' = 0. 

In summary, we have found in a strong coupling expan- 
sion to order (t/U) 4 that the energy gap in the insulating 
phase of the half-filled U-t-t' Hubbard model develops a 
(cos k x — cos k y ) 2 modulation at low temperatures. With- 
out the need to include hopping amplitudes beyond nnn 
an excellent fit is achieved for the single hole dispersion 
in S^CuC^C^. With the same dispersion fit parameters 
also the calculated gap modulation amplitude compares 
well with the ARPES data. It is natural to expect that 
these spectral features in the insulator carry over to the 
metallic, doped case; the consequences and the connec- 
tion to anisotropic pseudogap structures or even d-wave 
superconductivity remain yet to be understood. 
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